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Pressure of the Non-equilibrium O(N) $ 4 Theory in the Large N Limit 
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We calculate the off-equilibrium hydrostatic pressure for the O(N) 3> 4 theory to the leading order 
in 1/N. The present paper, the first of a series, concentrates on the calculation of pressure in 
the non-equilibrium but translationally invariant medium. The Jaynes-Gibbs principle of maximal 
entropy is used to introduce the relevant density matrix which is then directly implemented into 
dynamical equations through generalised Kubo-Martin-Schwinger (KMS) conditions. We show that 
in the large N limit use of Ward identities enables the pressure to be expressed in terms of two point 
Green's functions. These satisfy the Kadanoff-Baym equations which are exactly solvable, and we 
explicitly calculate the pressure for three illustrative choices of p. 
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I. INTRODUCTION 

In recent years significant progress has been made in un- 
derstanding the behaviour of QFT systems away from 
thermal equilibrium. Motivation for the study of such 
systems comes both from the early universe as well as 
from the quark-gluon plasma (deconfined phase quarks 
and gluons). Non-equilibrium effects are expected to be 
relevant in the relativistic heavy-ion collisions planed at 
RHIC and LHC in the near future 

One of the significant physical variables, in the con- 
text of non-equilibrium QFT, is pressure. Pressure, as 
an easily measurable parameter^], is expected to play an 
important role in a detection of phase transitions. This 
is usually ascribed to the fact that the pressure should 
exhibit a discontinuity in its derivative (s) when the local 
phase transition occurs. The aforementioned has found 
its vindication in solid state physics and fluid mechan- 
ics, and may play a crucial role, for instance, in various 
baryogenesis scenarios. 

It is well known that for systems in thermal equilib- 
rium, the pressure may be calculated via the partition 
function p|-|ll|. However, this procedure cannot be ex- 
tended to off-equilibrium systems where is no such thing 
as the grand-canonical potential. In this letter we con- 
sider an alternative definition of pressure, based on the 
energy-momentum tensor. This, so called, hydrostatic 



*In this connection we may mention the piezo resistive silicon 
pressure sensors used, for instance, in superfluid helium Q ^) 
or neutron (X-ray) diffraction technique used in solid state 
physics |,|. 



pressure is defined as the space-like trace of the energy- 
momentum tensor jl^ [l3| , and in equilibrium, it is iden- 
tified with the thermodynamical pressure via the virial 
theorem ||. There are several problems with the va- 
lidity of this identification on the quantum level, indeed 
gauge theories suffer from a conformal (trace) anomaly 
and require special care ||. However, we will avoid such 
difficulties by focusing on a scalar theory which is free of 
the mentioned complications ||[l3). The major advan- 
tage of defining pressure through the energy-momentum 
tensor stems from the fact that one may effortlessly ex- 
tend the hydrostatic pressure to non-equilibrium systems 
(for discussion see Ref. fu|)- 

The aim of this paper is to provide a systematic pre- 
scription for the calculation of the hydrostatic pressure in 
non-equilibrium media. This requires three concepts; the 
Jaynes-Gibbs principle of maximal entropy, the Dyson- 
Schwinger equations and the hydrostatic pressure. In 
order to keep our discussion as simple as possible we 
illustrate the whole procedure on a toy model system, 
namely the O(N) <E> 4 theory. The latter has advantage 
of being exactly solvable in the large N limit, provided 
that we deal with a translationally invariant medium. As 
a result the hydrostatic pressure may be expressed in a 
closed form. 

In order to provide meaningful results also for readers 
not entirely familiar with the Dyson-Schwinger equations 
and the Jaynes-Gibbs principle of maximal entropy, we 
briefly summarise in Sec. II the basic essentials. (More 
detailed discussion may be found in Refs. (l^Q.) As 
a byproduct we construct the generalised KMS condi- 
tions. Sec. Ill is devoted to the study of the (canonical) 
energy momentum tensor in the O(N) <I> 4 theory. If both 
the density matrix and the full Hamiltonian are invari- 
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ant under O(N) symmetry one obtains Ward's identities 
in a similar manner as in equilibrium. We show how 
these drastically simplify the expression for the pressure. 
In Sec. IV we concentrate our analysis on the large N 
limit. In this setting we derive a very simple expres- 
sion for the pressure - pressure depends only on two- 
point Green functions. Sec.V then forms the vital part 
of this paper. Owing to the fact that the infinite hier- 
archy of the Dyson-Schwinger equations is closed (basi- 
cally by chance) we obtain simple equations for two-point 
Green functions - Kadanoff-Baym equations. These are 
solved exactly for three illustrative density matrices p. 
We choose deliberately translationally invariant p's. The 
reason is twofold. Firstly, for a non-translationally in- 
variant medium one must use the improved energy mo- 
mentum tensor instead of the canonical one ||. This is 
rather involved and it will be subject of our future paper. 
Secondly, the Kadanoff-Baym equations turn out to be 
hyperbolic equations whose fundamental solution is well 
known. As a result we may evaluate, for the density ma- 
trices at hand, the hydrostatic pressure explicitly. The 
paper ends with a discussion. 



II. BASIC FORMALISM 



The key object of our interest is the energy-momentum 
tensor Q liV (x). A typical contribution to Q fiU (x) can be 
written as 



V^{x) 



(2.1) 



Here $ is a field in the Heisenberg picture and stands 
for a corresponding differential operator. Since 2? Mi $(a;) 
and 2?^ . $(x) generally do not commute for i j, one 
must prescribe the ordering in 8^,,. Our strategy is based 
on the observation that one can conveniently define such 
ordering via the non-local operator 



lim T*{0 Ml $(a:i)---2V*(a; n )} = 

lim ©(*a M )T{*(xi).. •*(«„)}, 



(2.2) 



where V(id^) is just a useful short-hand notation for 
D^D^ . . . 2? jLln , and T* is defined in such a way @ that 
for the product of fields without derivatives it is simply 
the T product, whilst for fields containing derivatives it 
is the usual T product with all the derivatives pulled 
out of the T-ordering symbol. It is clear that both T* 
and T coincide if all the arguments different, so 

(T* — T)(. . .) is an operator with a support at the contact 
points. The T* ordering is in general a very useful tool 
whenever one deals with composite operators. In the 
sequel we shall repeatedly use this fact. 



A. Off-equilibrium Dyson-Schwinger equations 

Let us now briefly outline the derivation of the Dyson- 
Schwinger equations for systems away from equilibrium. 
For simplicity we illustrate this on a single scalar field $. 

We start with the action where $ is linearly cou- 
pled to an external source J(x). Working with the fields 
in the Heisenberg picture, the operator equation of mo- 
tion can be written as 



SS 
8<§>(x] 



■[$ = $ J ] + J(x) = 0, 



(2.3) 



where the index J emphasises that $ is implicitly J- 
dependent. This dependence will be made explicit 
via a unitary transformation connecting $> J (governed 
by H — J$) with $ (governed by H). If J{x) 
is switched on at time t — ti we have $> J (x) = 
Tc{$(x) exp (i f c d 4 y J(y)$(y))}, where the close-time 
path C - the Keldysch-Schwinger path - runs along the 
real axis from ti to tf (tf is arbitrary, tf > U) a nd then 
back to ti. With this setting we can rewrite (2.3) as 



t; 



1 5S T> 



5$ 



J± 



xexp(i J d*y (J + (y)$ + (y) - J_(j/)$_(y)))j =0, (2.4) 

where, as usual || we have labelled the field (source) 
with the time argument on the upper branch of C as 
$+ (J+) and that with the time argument on the lower 
branch of C as $_ (J-)- Introducing the metric {(Jz) a i3 
(a 3 is the usual Pauli matrix and a, @ = { + ;—}) we can 
write - J_$_ = J a (ct 3 ) q/3 $f3 = J a Ma/J 

For the raised and lowered indices we simply read: = 
$+ and <!>_ = -$~ (similarly for J a ). Taking Tr(p. . .) 
with p = /?[<!>, 9$, . . .] being the density matrix, we obtain 
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SS 



Z[J] S$(x) 



$ Q (x) 



SJ a (x) 



Z[J] = -J a (x), (2.5) 



with Z[J] = Tr {p T c exp (i f c d A y J(y)$(y)) } being the 
generating functional of Green's functions. Employing 
the commutation relation: — i-gj-Z = Z {(f> a — igj-), we 



may recast (2.5) into more convenient form, namely 



J a (x) 



SS 
8®(x) 



+ z d 4 yG a (x iy ) (^ 



8<tP{y) 



1, (2.6) 



where the symbol 1 indicates the unit. As usual, the 
mean field, cf> a , is defined as the expectation value of 
the field operator: <j> a (x) = (& a (x)). Defining Z[J] = 
exp(iW[J]), two-point Green functions are then de fine d 
as G aP {x,y) = ~S 2 W/5J a {x)8J p {y). Setting J in @- 
to (i.e. physical condition) we obtain a first out of 
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infinite hierarchy of equ ation s for Green's functions. Suc- 
cessive J variations of (2.5)-(2.6) generate higher order 
equations in the hierarchy. The system of these equations 
is usually referred to as the Dyson-Schwinger equations. 

For the future reference is convenient to have the ex- 
pression for the effective action r[(/> c ]. This is connected 
with W[J] via the Legendre transform: 

T[ct> c ] = W[J] - f d A y J(y)My) ■ (2-7) 
Jc 

Following the previous reasonings, one can easily per- 
suade oneself that the expectation value of Q^ v reads 



+ i d 4 y G +/3 (x,y) (ct 3 ) 



Stiff (y) 



}1. (2. 



We have automatically used the sub-index '+' as the 
fields involved in Ml/ have, by definition, the time ar- 
gument on the upper branch of C. 



B. The Jaynes-Gibbs principle of maximal entropy 

Another important concept to be used is the Jaynes- 
Gibbs principle of maximal entropy. This prescription 
allows one to construct the density matrix in such a way 
that the informative content of p, subject to certain con- 
straints imposed by our theory/experiment, is the least 
jj^[l6|,[li|. The practical side of the principle rests on 
the maximalisation of the information (or Shannon) en- 
tropy, defined as S[p] — — Tr(p hip), subject to con- 
straints imposed by our knowledge of expectation values 
of certain operators Pj[$, 9$]. In contrast to thermal 
equilibrium, the Pj[. . .]'s need not to be the constants of 
motion (space-time dependences are allowed). 

In practice the Jaynes-Gibbs principle is used in the 
following way; given the constraints 

(P k [*,d*](x 1 ,...))=g h (x 1 ,...), k = l...n, (2.9) 

the maximum of the information entropy then determines 
the most plausible p as 



1 " 

P= zU7 exp(_ ^n rf4 ^ A*(si,...)flk[*,3*]), 

L ' fcJ i — 1 „• 

(2.10) 



fc=i j 



where the 'partition function' Z is 
Z[X k ] =Tr J cxpt-^IpV, \ k (x u ...)P k [<i>,d<i>}) 

{ fe=l 3 

In both previous expressions the time integration is either 
not present at all (so f k are specified only at the initial 
time ti), or is taken over the gathering interval (— r, £j). If 



instead, one has knowledge about the expectation values 
of /»[...] at some discrete times the corresponding 
integration must be replaced by a discrete summation. In 
the same manner as at equilibrium one may eliminate the 
Lagrange multipliers by solving the n simultaneous 
equations 



9k(xi, ■ 



5\nZ 



SX k (xi, .. .)' 



(2.11) 



The solutions of ( 2.11 ) may be formally written as 
_ S S G [gi, ...,g. 



Afc(xi, 



jn \ \max 



Sg k (xi, . . .) 



(2.12) 



Now, in order to reflect the density matrix ( 2.1C ) in the 
Dyson-Schwinger equations, we need to construct the 
corresponding boundary conditions. This may be done 
quite straightforwardly. Using the cyclicity of Tr(. . .) to- 
gether with the relation 

oo 

e A Be- A = V -C n , C = B,C n = [A, C^}, 
* — ' n\ 

n=0 

we can write the generalised KMS conditions for the n- 
point Green function as 

($(£!) . . . $(i„)) = ($(x 2 ) . . . $(a: n )$(a!i)) 

+ ^^{*(x 2 )...*(x n )C k ( Xl )), (2.13) 
fe=i 

where A = ln(p), B — $(xi), X\q — ti. So for the two- 
point Green functions we haveQ 

oo ^ 

G+-(x,y) = G-+(x,y) + V -Tr(p 4>{x)C n {y)). 

£ — ' TV. 

n—1 

In this paper and the paper to follow we aim to demon- 
strate that conditions (2.13) together with the causality 
condition are sufficient to determine Green's functions 
uniquely. 



III. THE O(N) $ 4 THEORY 
The O(N) <I> 4 theory is described by the bare Lagrangian 

JV . / N \ 2 

£ = 2((» a ) 2 -* 2 )^ £(<n 2 • 

a=l \a=l J 

(3.1) 



^In special cases when p — |0)(0| or p — e _/3H / ' Z{(3) the 
boundary conditions are the well known Feynman and KMS 
boundary conditions, respectively. 



3 



The corresponding canonical energy-momentu m t ensor 
is 6^ = J2a d fi $ a d"<f> a - g^C, and from Eqs.(fJ) and 



(2.8) its expectation value is 



9^ ($+(*)) = = *E 9 ^ G c + c + (x,y)\ x = y 



J2\^d ay G™ + (x,y)\ x -_ 



mlG% + {x,x) \ + 



8A 



c,d ^ ^ 



+ iG% + {x,x) ) ( (cf> d + (x)) 2 + iG d + d + (x,x) 



+ 2\^> c + {x)<j> d (x) + iG? + (x, x)j iGf + {x, x) 

+ terms with T 3 and T 4 . (3.2) 

Before proceeding further with our development , we want 
to show how one can significantly simplify Eq.(3.2) pro- 
vided that both the density matrix and the Hamiltonian 
are invariant under rotations in the TV-dimensional vector 
space of fields. This situation would occur if the system 
was initially prepared in such a way that no field $° was 
favoured over another. The fact that p is invariant under 
O(N) transformations means that 



U(e)p{<S>,d<S>,...]U- 1 (e)=p{<S>,d<S>,...} 



(3.3) 



where the fields <E> a transform under A-dimensional ro- 
tations: f/(e)$ a f/- 1 (e) = R- 1 (e) ab $ b = [exp(e l T l )] ab $ b , 
where R(e) is the rotation matrix in the iV-dimensional 
vector space, and the generators Ti are real and antisym- 
metric N x N matrices. It is obvious that the previous 
relation for <£ c can be satisfied for all times only if the 
full Hamiltonian, which governs the evolution of $ a , is 
also invariant against the A-dimensional rotations. 

Let us now consider the generating functional Z[J] cor- 
responding to the O(N) symmetric density matrix. Em- 
ploying the cyclic property of Tr(. . .) together with the 
infinitesimal transformation, SR(e) = 1 + e^X^, we obtain 
that the variation of Z must vanish. The latter implies 
the following (unrenormalised) Ward's identities: 



SJ b (y) 



= d*y J\y) 



5W[J] 
5J a (y) 



(3.4) 



Taking successive variations with respect to source J, 
we obtain the following results (see also (l6)): n-point 
Green functions with n odd vanish, while for n even (n = 
2k, k = 1, 2, . . .) one has 



<~>™, (Xl, . . . , X2k) 



ax---oi(2k) 



\^ TTa r (2fc) ( 
all pairings i<j 



,X2k) 



(3.5) 



where G^ 2k > is a universal 2fc-point Green function. Sim- 
ilar results fij| can be obtained for T^ 2 y a ^ k (...). 



Finally note that these results enable the expression for 
the expectation value of the energy momentum tensor to 
be simplified significantly to 

(e»"(x))=iNd>WG ++ (x,y)\ _„ 



N + 2 



\ g^{G ++ {x,x)f 



- l ^(d^ y G ++ (x,y)\ x=y 

- m 2 G ++ (x,x)^j + terms with T (4) . (3.6) 

In the rest of this paper we shall confine ourselves only 
to situations where both p and H are 0(N) invariant. 

IV. THE LARGE N LIMIT 



Let us now examine behaviour of (3J3) to the order 
0{1/N). For this purpose it is important to know how 
either G^ or r'™' behave in the TV — > oo limit. At 
T = or in equilibrium the Feynman diagrams are avail- 
able [L3| and the corresponding combinatorics can be 
worked out quite simply. On the other hand, the situa- 
tion in off-equilibrium cases is more difficult as we do not 
have at our disposal Wick's theorem. One may devise 
various diagrammatic approaches, e.g. kernel method 
Q , cumulant expansion M , correlation diagrams pi] , 
etc. Instead of relying on any graphical representation 
we show that for both equilibrium and off-equilibrium 
systems, the situation may be approached from far more 
general standpoint using only Ward's identities and prop- 
erties of T and W. 

In order to find the leading behaviour at large N it is 
presu mab ly the easiest to consider the Legendre trans- 
form ( |2.7| ) . The explicit N dependence may be obtained 
by setting 4> c = </>, which implies J c = J for all the group 
indices. Eq.(2.7) then indicates that both T[<f>] and W[J] 
are of order N. If we expand T[<j> a } in terms of <j> c taking 
into account Ward's identities we get 

m = r[0] + In [ d A x d A y T (2 \x,y) <f>(x)<f>(y) 



+ ^n 2 [ AdV 4 2Ar (4) (i,!/,M) 

4! Jc 



x (j)(x)(/)(y)(j)(z)(j)(q) 



(4.1) 



Since the LHS of Q is of order AT, must be of 
order TV , T^ 4 ) of order iV _1 , and, in general, r^ 2 ") must 
be of order iV 1_ ". This suggests that in the expression 
for the energy-momentum tensor ( |3.6[ ), terms containing 
r^ 4 ) can be ignored. The above argument can be repeated 
in a similar way for W. 

Hence, collecting our results together, the expectation 
value of the energy-momentum tensor to leading order in 
N may be written as 
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= iNd%d» G ++ (x,y) 

.N 
~2 



i-g^id^d" G ++ [x,y)\ -ml G ++ (x,x) 



(4.2) 



This result is surprisingly simple: the expectation value 
of the energy-momentum tensor, and thus the hydro- 
static pressure, is expressed purely in terms of two- 
point Green's functions. The latter can be calculated 
through the Dyson-Schwinger equations ( |2.6| ). Further- 
more, these equations have a very simple form provided 
that both the large TV limit and Ward's identities are ap- 
plied. If we perform a variation of ( [2.6| ) with respect to 
Jp(y) we obtain, to order 0(1/N), the following Dyson- 
Schwinger equations for two-point Green functions: 



□ 



+ ™ 2 + iX ° r 

' m "I ^aa 

-8(y - x)(<T 3 ) a p . 



(x,x) G a p(x,y) 



(4.3) 



These dynamical equations for two-point Green functions 
are better known as the Kadanoff-Baym equations . 

Let us mention one more point. The generalised KMS 
conditions for G± T are significantl y sim ple in the large 
N limit. This is because in sum ( 2.13] ) only terms of 
order O(N ) contribute. This implies that only quadratic 
operators Pi [$, 9$] in the density matrix are relevant. As 
a result, the Jaynes-Gibbs principle naturally provides a 
vindication of the popular Gaussian Ansatz 



V. OUT-OF-EQUILIBRIUM PRESSURE 

The objective of the present section is to show how the 
outlined mathematical machinery works in the case of 
the hydrostatic pressure. In order to gain some insight 
we start with rather pedagogical, but by no means triv- 
ial examples; translationally invariant, non-equilibrium 
density matrices. We consider the more difficult case of 
translationally non-invariant density matrices in a future 
paper. 



A. Equilibrium 

As an important special case we can choose the con- 
straints 



gk = constant, 



(5.1) 



where t% is arbitrary. Eq. fl5.ip then implies that iVs 
are integrals of motion. Since in the finite volume sys- 
tems the spatial translational invariance is destroyed, the 
only integral of motion (apart from conserved charges) is 
the Hamiltonian. Thus the system is in thermal equilib- 
rium and the laws of thermodynamics pi prescribe that 



g = L dT'Cv(T') (Cv is the heat cap acity at constant 
volume V and T is temperature). Eq. (2.12) determines 
the Lagrange multiplier; A = 1/T = (3. The density ma- 
trix maximising the Sq is then the density matrix of the 
canonical ensemble: p = — -• Due to the transla- 
tional invariance of p, the Kadanoff-Baym equations read 



(□ a + m 2 r (T)) G a(3 (x -y) = -S(x - y)(a 3 ) a(3 , (5.2) 

where the temperature-dependent renormalised mass is 
(see, for example [HI El); m 2 r {T) = m 2 Q + ^ G aa (0). 



(5.3) 



The corresponding KMS boundary condition is 
G+_ (x, U ; y , 0) = G_+ (x, U - i(3; y , 0) . 



Because m r (T) is a spatial constant, a Fourier transform 
solves equations (p^). The solutions of (5.2) subject to 
condition ( |5.3| ) arc then the resumed propagators in the 
Keldysh-Schwinger formalism 

tG ±T (k) = 2Tr{0(Tko) + f(\k o \)}5(k 2 -m 2 r ) , (5.4) 

with /(x) = (exp(/3x) — 1) _1 being the Bose-Einstein 
distribution. 

Now, the total hydrostatic pressure in d dimensions is 
classically defined as ||,[l|,[L3| 



V{T) 



(d-1) 



Because is a composite operator, a special renor- 
malisation is required |,[|||,|2|. As we have shown 
in Jl5|, for translationally invariant systems the renor- 
malised pressure coincides with the, so called, thermal 
interaction pressure P. The latter reads 



P(T) = V(T) - P(0) 

= -j^ T) {m-m i M} 



(5.5) 



Result (5.5) deserves two comments. Firstly, the energy 
momentum tensor 0^ need not to be the canonical one 
(however, the canonical one is usually the simplest one). 
This is because energy momentum tensors are mutually 
connected via Pauli's transformation, i.e. 



0^ = 



W + d x X x ^ 



(5.6) 



and owing to the space-time translational invariance 

(©>") = (6^> + d x {X x ^) = (0^). (5.7) 

(The analogous identity is naturally true at T = 0.) 
As a second point we should mention that prescrip- 
tion d5.5| ) retains its validity for non-equilibrium media 
as well. This is because the short-distance behaviour 
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of G++(x, y), which is responsible for the singular be- 
haviour of M ", comes fro m the particular solution of the 
Kadanoff-Baym equation (5.2). The latter can be cho- 
sen to be completely independent of the boundary con- 
ditions (actually it is useful to chose the Feynman causal 
solution). On the other hand, the homogeneous solution 
of (j5~2j), which is regular at \x — y\ — > 0, reflects all the 
boundary conditions. One may see then that the UV sin- 
gularities which trouble 8 Ml/ may be treated in the same 
way as at the T = 0. Incidentally, the former is an exten- 
sion of the well known fact that in order to renormalise 
a finite temperature QFT, it suffices to renormalise it at 
T = 0. 

Inserting the solution ( b-jf) i nto the expression for the 
energy-momentum tensor (|4. 2J ) we arrive at the thermal 
interaction pressure per particle (see also |Ia,p§|) 



P{T) = 



T 2 m 2 r {T) Tm z r {T) 



90 



24 



12tt 



rp4 



T 3 m r {T) T 2 m 2 r (T) 



24tt 



16tt 2 



+ O m 4 (T) In 



m r (T) 
T4n 



(5.8) 



where the renormalised coupling constant A r comes from 
the T = renormalisation prescription: r^2_ >bb (s = 0) = 
— \ r /N (s is the usual Mandelstam variable). A direct 
calculation of I^ 4 ) in the large N limit |ii],[k|, gives the 
relation for the renormalised coupling: 



A 



2^(4^72 



r[2-d/2] K) 



2\d/2-2 



(5.9) 



To regularise the theory we have used the usual MS 
scheme with the dimensional regularisation (d 4 ). The 
high temperature expansion of the pressure ( p.8[ ) to all 
orders can be found in |Q where the calculations, how- 
ever, were approached from a different standpoint. 



B. Off-equilibrium I 

The next question to be addressed is how the above cal- 
culations are modified in the non-equilibrium case. To 
see that let us choose the following constraint 



S (k) = (W(k))| t4 = (W(k))| t4 . 



(5.10) 



Here ti is arbitrary, and function g(k) is specified by the- 
ory/experiment. The 7i(k) is a quadratic operator ful- 
filling the condition (H.)\n^oo = (W). As g(k) is finite, 
both Ti and Ti must be renormalised (i.e. we must sub- 
tract the zero temperature part). Obviously 

H(k) =« fc ot(k)a(k), 

with ujk = Vk 2 + M 2 and 



M 2 = ml 



1 x f ddk ' r 



(5.11) 



In the large N limit the corresponding density matrix 
Q2.10D reads 



exp 



d 3 k 



{2ir) 3 2iu k 



/J(k)W(k) 



(5.12) 



with /3(k)/2(27r) 3 a;fe be ing th e Lagrange multiplier to be 
determined. Using Eq.(2.11), we find that 



ff(k) = 



V 



(2tt) 3 e/W"* - 1 ' 



(5.13) 



where V deno tes the volume of the system. Clearly, ex- 
pression (5.13) can be interpreted as the density of energy 
per mode. The fact that /3(k) is not constant indicates 
that different modes have different 'temperatures'. 

The Kadan off- Baym equations coincide in this case 
with those in (5.2) provided m r (T) — > A4. The boundary 
condi tion can be worked out simply using the prescrip- 
tion fl2.13|). This gives 



G+_(fc) =e- p ^ ko G- + {k). 



(5.14) 



The fundamental solution of the Kadanoff-Baym equa- 
tion is 



G af} {k) 



(03)a 
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k 2 -M 2 



2niS(k 2 -M 2 )f a p(k), 



tap = e(0"3) a /3 . 



(5.15) 



Let u s mention one more point. The boundary condition 
([5.14 ) is not by itself sufficient to determine / Q( 3's. (This 
fact is often overlooked even in the equilibrium QFT.) It 
is actually necessary to substitute this condition with an 
additional condition, namely the condition of causality. 
The causality condition, i.e. vanishing of the commuta- 
tor [$(x), $(y)] for (x — y) 2 < 0, importantly restricts 
the class of possible / a/ a's. To see this, let us look at the 
Pauli- Jordan function ( The latter is the 
homogeneous solution of the Kadanoff-Baym equation 
with the initial-time value data: ( [$(x), ®{y)])\x =y — 
(i.e. causality condition) and do{ [&(x), ®(y)])\x =y = 
— £ 3 (x — y). Thus the Pauli- Jordan function is uniquely 
determined and its Fourier transform reads 

F.T.«[*(x), *(»)]» = -i2irS{k 2 - M 2 ) e(k ). (5.16) 



Relation (5.16) immediately implies that 



f + .(k) = e(-k ) + f(k) 

/_+(*:) = <?(*<>) + /(*), 



(5.17) 



with / being, so far arbitrary, and for both /_| and / |_ 

identical, functio n. T he / is then fixed via the generalised 
KMS condition (5.14). Similarly, the causality condition 
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specifies that G++(fe)-GL_(fc) =PP{l/(k 2 -M 2 )} (the 
symbol 'PP' denotes the principal part). By inspection 
of the definition of G a $ one may easily realise that 



G-\--\- -\- G — — G-\ — -\- G— 

G+-(k) = -{G+-(-k)y 
G++(k) = -(G—(-k)y 



(5.18) 



From these relations follows that / ++ = / = / and 

f(k) = (f(-k))*. The / is the same as in ( pT7j ). 



Since the divergences in G a p come from the first term 



in ( 5.15| ) (i.e. from the particular solution), we can shift 
the corresponding (zero temperature) poles at d — 4 to 
the bare mass. In this case we can write 



M 



5m 



(5.19) 



where m r is the renormalised mass in the vacuum and 
dm is the mass shift due to an interaction with the non- 
equili brium med ium. Inserting the 'H — h' components of 
(5.15) into ( |5.11 ), we obtain 



M 2 = m 2 + ^Ao 



M 



d-4 



-r [1 - d/2] + N(M 2 ) 



(5.20) 



where 



N(M 2 ) = 



d d k 



(27T) 



2 2nS{k 2 -M 2 )f ++ {k). (5.21) 



In an equilibrium system / ++ would be the Bose-Einstein 
distribution. Note that because (5.21) corresponds to a 
homogeneous solution of the Kadanoff-Baym equation at 
\x — y\ = it is automatically finite. Thus for a trans- 
lationally invariant medium (both equilibrium and non- 
equilibrium) /_|__|_ must act as a regulator in the UV re- 
gion. 

Let us consider the expression for the expectation value 
of 0^ given in ( |4.2| ). It is a matter of a few lines to show 
that 



(9H 



iN 
,N 



d d k 
{2n) d 



k»k v [C++ (A) - G{k)] 
d d k 



(2tt) c 



[G++(k) + G(k)} , (5.22) 



with G be ing t he T — causal Green function. Us- 
ing (|5.9|) , (5.11) and (5.15) one may directly check that 
Eq. fl5.22[ ) does not contain UV singularities when the 
limit d — > 4 is t aken. This verifies our introductory re- 
mark, that ( |5.5|) is finite even in non-equilibrium con- 
text. From the generalised KMS condition (5.14) and 



from ( 5.1 5| ) we obtain that / ++ is 



/++(*) 



1 



(5.23) 



So far the results obtained were completely general and 
valid for any translationally invariant system. Let us now 



consider a system in which g(k) = X%f$ exp (— ujk/c). As 
we shall see, this condition corresponds to systems where 
the lowest frequency modes depart from strict equilib- 
rium, whilst the high energy ones obey standard Bose- 
Einstein statistics. This behaviour is typical of many 
off-equilibrium systems, eg. ionised atmosphere , laser 
stimulated plasma, hot fusion etc. The a is usually 
referred to as the ionization half-width of the energeti- 
cal spectrum. From ( 5.13 ) we can determine /3(k) as a 
function of the physical parameter a: 



1 1 00 
/3(k) = i + ^^ 
cr UJ k 



(-1) 



n+l 



(5.24) 



To proceed, some remarks on the interpretation of /5(k) 
are necessary. Firstly, Eq.(5.24) implies that for a suffi- 
ciently large LOk (y>k S> cr) the function /3(k) is approxi- 
mately constant, and equals l/cr. Thus at high energies 
the distribution / ++ approaches the Bose-Einstein dis- 
tribution with the global temperature T ~ cr. In other 
words, only the soft modes were sensitive to processes 
which created the non-equilibrium situation. The inter- 
pretation of cr as an equilibrium temperature, however, 
fails whenever crsiwj.. Instead of cr one may alternatively 
work with the expectation value of /3(k), i.e. 



Jd 3 k(3(k) 



i(n+l) 



1 , EZ^-^QK^n+^M/a) 



MK 2 (M/a) 



(5.25) 



where K n is the Bessel function of imaginary argument 
of order n. Because the system is for the significant part 
of the energetical spectrum in thermal equilibrium, l/(/3) 
approximates the corresponding equilibrium temp erature 
to high precision. An interesting feature of ( 5.25| ) is that 



it is quite insensitive to the actual value of M . Depen- 
dence of (/?) on M. for a sample value cr = lOOMeV is 
depicted in FIG.|l|, 



(P) [MeV] 




. 0104 



0.0102 



M [MeV] 



200 400 600 800 1000 

FIG. 1. A plot of (J3) vs. M atu= 100 MeV. 
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An important observation is that the asymptotic be- 
haviour of (P) at (7 — > oo goes like (f3) ~ a/c, where 

Using 
ticle 



) and (5.22) we get for the pressure per par- 



P = 



1 

2^2 



M 2 a 2 K 2 (M/a) 



5m 



-MaKi(M/a) 



-Vo, 
(5.26) 



where 



Vo 



1 



384tt 2 
1 

y 

64tt 2 ' 



(Sm 2 )m 



1 





(5m 2 \ 




m£ / 


<5m 2 N 


1 In ^1 ^ 


m 2 , 





(5m 2 



(5.27) 



Note that "Po comes from the UV divergent parts of 
( |5.22 ). Whilst the separate contributions are UV diver- 



gent, they cancel between themselves leaving behind the 
finite Vo- As we already emphasised, the divergences 
come from the particular solutions of the Kadanoff-Baym 
equations. Because the former do not directly reflect the 
boundary conditions, the form of Vo must be identical 
for any translationally invariant medium. The non-trivial 
information about the non-equilibrium pressure is then 
encoded in terms in the brace [. . .] in ( 5.26| ). 

For a sufficiently large o the leading behaviour of dm 2 
may be easily evalu a ted. To do this, let us first assem- 
ble d^|) and d5.19| )- (|5.2l]) together. This gives us the 
(renormalised) gap equation 



5m 2 



^^-TrJ^iM^-rnt^+NiM 2 ) 

\f' \ ~ , n _ o > 1 



^-{*£(m 2 r ,5m 2 ) + — 2 -aMKi[M/a] }, (5.28) 



with 



32^ 



(m 2 + (5m 2 )ln 1 + 



— Smt 



Setting x = 8m 2 jm 2 and s = a/m r we obtain the fol- 
lowing transcendental equation for x 



1 



32tt 2 x 



- [(1 + ar)m(l + x) - x 



+ 8sVTT^ K^l + x/s] ] } 



(5.29) 



The graphical representation of Eq. ( 5.29[ ) is depicted in 
FIG.|. From this one may read off that for large x (x > 
50) there exists a critical value of A r above which the 
gap equation has no solution. (The plateau is actually 
bent downward with a very gentle slope.) FIG.|2|b clearly 
shows that if A r 1/s < 1 then x -C 1. Using the 
asymptotic behaviour of Ki[z] for z — ► (K\ ~ (z) -1 ) 



we arrive at more precise estimate of A r for which i« 1; 
namely A r » 1/s 2 = m 2 /cr 2 . This estimate is very helpful 
for the asymptotic expansion of <5m 2 . For a sufficiently 
high a we may write 



5m 



\ r cr 



0(Mln(M/a)). 



(5.30) 



Inserting this result into ( ^3.26| ) and using the series repre- 
sentation of both K 2 and Ki, keeping the leading terms, 
we get 



P(a) 



(x 4 a 2 M 2 X r ( a 2 M 2 a 4 3\ 
7T 2 2tt 2 + 8 I 64tt 4 4tt 4 J 



+ 0(M 2 HM/(i);X 2 r ) 



(5.31) 




FIG. 2. A plot of the Eq. ( 5.2i ): a) the general shape, b) a 
small x behaviour. 



It is interesting to note that I/71" 2 w 7r 2 /97 which is 
almost the Stefan-Boltzmann constant. This once more 
vindicates our interpretation of a as a "temperature" . A 
plot of the pressure as a function of a is depicted in FIGj|. 
It is due to the low frequence modes contribution to the 
pressure that P(cr) < P(T). This is a direct result of our 
choice of g(k), namely that a cannot be interpreted as 
temperature for the low frequence modes (i.e. LOk < o~). 
The smaller a is the less important contribution from 
non-equilibrium soft modes and so the smaller difference 
between both pressu res. 

The result (5.31) can be alternatively rewritten in 
terms of (8) . Using, for instance, the Pade approxi- 
mation pel for (/3), we arrive at 
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P({/3)) = 0.0681122 (/3)- 4 - 0.0415368 (/3)~ 2 M 2 

+ \ r (-0.000647 (P)- 4 + 0.0000164 ((3)- 2 M 2 ) 
+ 0(M 2 \n(M(f3));\ 2 r ) . (5.32) 

The coefficient 0.0681122 w tt 2 /145 is 1.6 times smaller 
than the required value for the Stefan-Boltzmann con- 
stant, so the parameter ((3) is a slightly worse approxima- 
tion of the equilibrium temperature than a. In practice, 
however, it is a matter of taste and/or a particular con- 
text which of the above two descriptions is more useful. 



4 -10 
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FIG. 3. A plot of pressure as a function of T, a for 
m r ~ WOMeV. The gray line corresponds to equil ibrium pres- 
sure, the black line corresponds to pressure (5.31). 
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FIG. 4. A plot showing the difference of equilibrium and 
non- equilibrium pressures for m r — lOOMeV. 



'slightly' from equilibrium, i.e. when g(k) in ( 5.1 0| ) de- 
viates a little from the equilibrium density of energy per 
mode. In this case the constraint ( 5 . 1 C| ) reads 



.9(k) 



V 



(2tt) 3 e^h - 1 



Sg(k); 5 3 (k)« 5 (k), (5.33) 



with /3q = 1/Tq being an inverse of the equilibrium tem- 
perature. As a special example of ( |5.33 ) we choose 



.9(k) 



V 



(2tt) 3 



a(k,A))wl. (5.34) 



The inverse mode "temperature" /3(k) is then (3o — 
ln(a(k))/u>k- So ln(cc) measures (in units of u) k ) the devi- 
ation of the mode temperature from the equilibrium one. 
The particular value of a(k, /3 ) depends on the actual 
way in which the system is prepared. To avoid unnec- 
essary technical complications, we select a(k, /?o) to be 
a momentum-space constant (generalisation is, however, 
straightforward). This choice represents the change in 
the mode temperature which is now inversely propor- 
tional to ui k \ the deviation is bigger for soft modes and 
is rapidly suppressed for higher modes. Obviously, To 
becomes the global temper atur e if u> k 3> ln(ct). The gen- 
eralised KMS conditions ( 5.14 ) together with solutions 
( |5.15 ) determine /++ as 



/++(*) 



a 



(5.35) 



Eq.( |5.35|) is a reminiscent of the, so called, Juttner dis- 
trjbution^with fugacity a [pi|,p^|. Now, using fl5.5| ) and 
fl5.22fl we get for the pressure per particle 



KA2 00 n 

. n=1 n 
MSm 2 ^ a n 
4/3 ^ n 

where Sm 2 satisfies the gap equation 
dm = — \±j(m r ,dm ) 



(5.36) 



27r 2 J u k e^°" fc — a t 

(5.37) 



If we set, as before, x = Sm 2 /m 2 and s = l/pom r we get 
the transcendental equation for x 



C. Off-equilibrium II 

As was noted above, it is the specific form of the con- 
straint g(k) which prescribes the behaviour of /3(k). 
Let us now turn our attention to systems which depart 



'it should be mentioned that this similarity is rather su- 
perficial. The Jiitner distribution describes systems which 
are in thermal but not chemical equilibrium. (As we do not 
have a chemical potential, chemical equilibrium is ill defined 
concept.) The fugacity then parametrises the deviation from 
chemical equilibrium. 
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a 



T [GeV] 



+ 8(1+ x) 



dz \fi 



3 z y'x+l/ s 



(5.38) 



The corresponding numerical analysis of ( 5.3§| ) reveals 
that for x < 1 , 1/s < 1. So at i « 1, A,, w <5m 2 /T 2 . 
This estimate is important for the asymptotic expansion 
of 5m 2 . Howe ver, in order to carry out the asymptotic 
expansion of ( 5.37 ) (and consequently ( 5.36 )) we need 
to cope first with the sum X^^Li ot n Kk(nM.(3o) /n k (k 
1.2 (also called Braden's function). Expansion of 
Kk{. ■ .) yields a double series which is very slowly con- 
vergent, and so it does not allow one to easily grasp the 
l eadin g beha viour in To- ln this case it is useful to resume 
( 5.36 ) - J5I37 ) by virtue of the Mellin summation technique 
U . (It is well known |)], |l3|, |3J| that at equilibrium this 
resummation leads to a rapid convergence for high tem- 
peratures.) As a result, for a sufficiently large T we get 



(5m 2 



X r r 2 \ r T M 
24 



V-iry 
M 



I ( 1-^)V2 



/ I 1 - ln ( ])-(!- r 2 ) 1/2 arcsin(r) 



O(lnTo) , 



(5.39) 



where we have set r = hx (a)Tg /M.. The corresponding 
expansion of the pressure ( |5.36| ) reads 



P 



90 



^C(3) s r 



- 5 (l-r 
|(1 -<> 



24 

-It- -r' 



r 



2 

3 , 3 



3 
3(1 



(1 - 2, 2 ) - 
1 - 2 ln 



M 3 T 
4tt 2 
M 
2T 



2\ 2 



2 ) 



(5m 2 ) 2 

2A r 



■O(lnTo) 



arcsin(r) 

(5.40) 



where £(3) « 1.202. N ote that for a = 1 we regain the 
equilibrium expansion (|5.8| ). The corresponding plot of 
P as a function of T and a is depicted in FIG.H. 




FIG. 5. Behaviour of the pressure (5. St ) as a function of 
a and T at m r = lOOMeV. 
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In passing it may be mentioned that the expansion ( 5.40| ) 
is mathematically justifiable only for a e V A >'/ 24r « l. 



VI. SUMMARY AND CONCLUSIONS 

In order to get a workable recipe for non-equilibrium 
pressure calculations we have combined three indepen- 
dent concepts: the Jaynes-Gibbs principle of maximal 
entropy, the Schwinger-Dyson equations, and the hydro- 
static pressure. The basic steps are as follows. 

To find quantitative results for pressure one needs to 
know the explicit form of the Green's functions involved. 
These may be find if we solve the Dyson-Schwinger equa- 
tions. The corresponding solutions are unique provided 
we specify the density matrix p (the construction of p is 
one of the cornerstones of our approach, and we tackled 
this problem using the Jaynes-Gibbs principle of maximal 
entropy). With p at our disposal we showed how to for- 
mulate the generalised Kubo-Martin-Schwinger (KMS) 
boundary conditions. 

To show how the outlined method works we have il- 
lustrated the whole procedure on an exactly solvable 
model, namely 0(N) <i> 4 theory in the large N limit. 
This model is sufficiently simple yet complex enough to 
serve as an illustration of basic characteristics of the pre- 
sented method in contrast to other ones in use. In or- 
der to find the constraint conditions we have considered 
two gedanken experiments in which the system in ques- 
tion was prepared in such a way that only low frequence 
modes departed from the strict equilibrium behaviour. 
Such processes can be found, for example, in ionised at- 
mosphere, in laser stimulated plasma or in hot fusion 
Q. In both alluded gedanken experiments we were 
able to work out the hydrostatic pressure exactly. Fur- 
thermore, after identifying a "temperature" parameter 
(virtually temperature of high modes) we carried out the 
corresponding high-temperature expansions. 

As it was discussed, one of the main advantages of 
the Jaynes-Gibbs construction is that one starts with the 
(physical) constraints (i.e. parameters which are really 
controlled and measured in experiments). These con- 
straints directly determine the density matrix and hence 
the generalised KMS conditions for the Dyson-Schwinger 
equations. This contrasts with the usual approaches 
where the density matrix is treated as the primary object. 
In these cases it is necessary to solve the von Neumann- 
Liouville equation. This is usually circumvent using ei- 
ther variational methods |2^, |l| with several trial p's 
or reformulating the problem in terms of the quantum 
transport equations for Wigner's functions It is, 

however, well known that the inclusion of constraints 
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into transport equations is a very delicate and rather 
complicated task (the same is basically true about the 
variational methods) 
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